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ABSTRACT 


We^show  that  many  of  the  ^false  alarms'^  which  require  analyst  intervention 
in  the  preparation  of  the  LASA  event  summary  were  due  to  local  or  regional 
events.  Analyses  showed  that  these  false  alarms  occur  predominantly  on  weekdays 
during  local  working  hours,  suggesting  that  the  seismic  events  are  of  man-made 
origin.  The  false  alarm  ra*-e  decreases  on  weekends  and  holidays,  and  LASA 
reports  more  teleseismic  ev  'ts. 

To  reduce  the  number  of  false  alarms  it  is  necessary  to  steer  detection 
beams  to  local  areas.  By  detecting  local  events  on  these  beams  and  by  using 
a higher  S/N  threshold  in  processing,  these  signals,  we  can  ef-f actively _x;educe 
the  number  of  false  alarms'jj^from  the  original  57%  to  41%. 

A new  beam  set  has  been  developed  and  deployed  which  concentrates  tele- 
seismic  beams  in  high  seismicity  areas  instead  of  spacing  them  equidistantly 
apart.  This  arrangement  reduced  the  average  detection  errors  from  200  km  to 
50  km,  there  is  also  some  indication  of  a lowered  detection  threshold  on  the 
order  of  0. 1^0.1  magnitude  units. 
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INTRODUCTION 


The  on-line  data  processing  of  LASA  short-period  data  at  the  Seismic 
Data  Analysis  Center  (SDAC)  was  operational  from  January  1971  to  June  1975. 
During  this  period  of  operation,  the  main  purpose  was  to  produce  promptly  and 
routinely  an  SDAC/LASA  daily  event  summary.  The  operation  of  SDAC  was  tempo- 
rarily suspended  in  July  1975  in  order  to  adjust  to  the  instrument  and  data 
reconfigurations  at  LASA,  and  to  implement  modifications  which  would  enable 
SDAC  to  accept  and  process  seismic  data  from  additional  stations. 

Data  processing  is  performed  in  two  parts  by  utilizing  two  computers. 

The  first  part  is  the  Detection  Processor  (DP)  which  performs  on-line  signal 
detection  by  forming  and  applying  detection  algorithms  to  a number  of 
surveillance  beams.  The  second  part  is  the  Event  Processor  (EP)  which  selec- 
tively processes  detected  signals  and  extracts  event  parameters  such  as  event 
location,  origin  time  and  magnitude. 

The  term  "false  alarm"  has  been  historically  used  to  depict  noise  detec- 
tion! . False  alarm  in  this  sense  means  detections  on  fluctuations  of  noise, 
and  -an  include  instrumental  noise  such  as  transmission  errors.  Efforts  to 
reduce  this  type  of  false  alarm  have  been  made  in  the  past.  Lacoss  (1972) 
argued  that  since  false  alarms  are  detections  of  noise  fluctuations,  both 
noise  level  and  noise  variance  can  affect  the  rate  of  false  alairms.  Studies 
of  noise  variance  and  its  effect  on  signal  detectability  have  been  performed 
at  NORSAR  (Bungum  and  Husebye,  1974;  Bungum  and  Ringdal,  1974;  and  Steinert, 
Husebye  and  Gjoystdal,  1975).  In  particular,  Steinert  showed  that  the  noise 


Lacoss,  R,  T.,  1972,  Variation  of  false  alarm  rates  at  NORSAR;  Semiannual 
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Steinert,  0.,  E.  S,  Husebye,  and  H.  Gjoystdal,  1975,  Noise  variance  fluctuations 
and  earthquake  detectability;  Geophys.  Jf.  R.  Astr.  Soc.,  v.  41,  p.  289-302. 


stability,  which  is  the  measure  of  the  ratio  of  noise  average  to  noise  variance, 
is  Indeed  the  most  effective  Indicator  of  .he  false  alarms.  The  NORSAR  DP  is 
currently  operating  with  varying  detection  thresholds  based  on  the  noise  level 
and  noise  stability.  Chang  (1974)  argued  that  if  the  noise  fluctuations  can 
be  considered  as  random  occurrences,  then  the  requirement  for  a number  of 
consecutive  threshold  crossings  would  effectively  reduce  such  false  alarm 
detections.  In  comparing  LASA  and  NORSAR  detection  algorithms,  Chang  found 
that  the  temporal  requirement  of  consecutive  threshold  crossings  was  set  to 
one  at  NORSAR,  but  LASA  DP  required  three  consecutive  crossings.  In  any  case, 
study  of  past  operations  at  LASA  confirms  thr*-  noise  detections  were  indeed 
rare  at  LASA  and  did  not  pose  a significant  problem  to  the  DP  and  EP  operations. 

Chang's  comparison  of  LASA  and  NORSAR  short  ^eriod  array  performances 
shows  that  in  both  arrays  about  12%  of  the  DP  defections  are  ultimately 
published  as  events  in  dally  event  summaries.  This  of  course  does  not  mean 
all  DP  detections  are  processed  bv  EP.  Of  all  DP  detections  approximately 
one  third  were  processed  by  EP.  An  initial  reduction  of  DP  detections  is  made 
in  EP  by  a higher  threshold  setting  and  by  a grouping  algorithm.  Our  atten  ^'n 
is  drawn  to  the  fact  that  of  all  signals  completely  processed  by  EP,  only  ^ c 
half  of  them  were  confirmed  by  an  analyst  and  published  in  the  daily  summr 
In  the  early  evaluation  of  the  LASA/SAAC  system.  Dean  (1972)  reported  that 
only  37.3%  of  EP  processed  signals  were  reported  on  the  LASA  Dally  Summary. 

If  half  of  the  signals  processed  by  EP  are  false  (EP  false  alarms),  the  nature 
of  these  signals  should  be  investigated.  Since  it  seems  clear  that  DP  detec- 
tions are  relatively  free  of  noise  detections,  we  conclude  that  EP  false 
alarms  are  seismic  signals  that  cause  difficulty  in  the  processing  and  produc- 
tion of  the  event  bulletin. 

In  the  future  SDAC  Network  Event  Processor,  carefully  selected  signals 
of  one  station  will  be  associated  with  detections  of  other  stations.  The 


Chang,  A.  C.,  1974,  A comparison  of  the  LASA-NORSAR  short-period  arrays: 
SDAC-TR-74-5,  Teledyne  Geotech,  Alexandria,  Virginia. 

Dean,  W.  C.,  1972,  A geophysical  evaluation  of  the  short-period  LASA/SAAC 
system:  SAAC  Technical  Report  No.  5,  Teledyne  Geo tech,  Alexandria, 

Virginia. 


carefully  selected  signals  are  those  detections  processed  with  some  type  of 
process  which  will  be  in  fact  identical  to  the  EP  process.  Therefore  if 
these  siganls  contain  many  false  alarms,  it  will  be  difficult  to  obtain  good 
results.  In  this  study  false  alarms  are  classified  into  several  categories 
and  analyzed  to  show  the  rate  of  occurrence  of  each  type.  Discussions  of 
how  false  alarms  cause  detections  and  possible  methods  to  reduce  them  are 
tested  with  the  on-line  detection  processor. 


CAUSE  OF  FALSE  ALARMS  IN  THE  EVENT  PROCESSING  SYSTEM 


i 
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The  detection  threshold  of  the  LASA  DP  processor  has  been  set  to  10  dB. 

With  this  threshold,  there  are  approximately  300  detections  per  day.  Not  all 
of  these  detections  are  processed  by  EP.  A reduction  in  the  number  of  detections 
in  EP  depends  mainly  on  two  conditions:  thresholding  and  grouping.  The  first 

condition  simply  raises  the  processing  threshold  to  14  dB,  thus  eliminating 
approximately  two  thirds  of  the  detections.  The  remaining  100  detections  are 
then  screened  by  the  grouping  algorithm.  The  grouping  algorithm  checks  each 
detection  and  searches  for  consecutive  detections  that  are  restricted  in  area 
and  in  a specified  time  window.  Only  the  first  detection  in  the  grouped  detec~ 
tions  is  processed.  This  algorithm  further  reduces  by  30%  the  number  of 
detections  reaching  EP.  The  detections  screened  through  these  two  conditions 
are  processed  by  EP  on  a routine  daily  basis.  There  are  approximately  70  of 
such  events  per  day  that  reach  EP  for  final  processing. 

The  task  of  the  analysts  is  to  examine  results  of  the  automatically 
processed  EP  events  by  displaying  waveforms  and  seismic  parameters  on  the 
Experimental  Operations  Console  (EOC).  The  analyst  can  confirm,  adjust,  submit 
for  reprocessing,  or  reject  the  processed  event.  The  final  results  are 
published  in  the  SUMMARY  OF  SDAC/LASA  VELOCITY-BEAM  LOCATIONS,  which  contains 
an  average  of  30  events  per  day.  The  analyst  therefore  rejects  more  than  50% 
of  the  events  that  reach  EP  in  daily  operations.  These  rejected  events,  the 
EP  false  alarms,  are  the  main  interest  of  the  current  analysis. 

A data  period  of  37  days  from  24  June  to  30  July,  1974,  was  selected  and 
all  EP  processed  events  were  studied  and  grouped  into  seven  categories.  These 
seven  categories  are;  (1)  Identified  events,  (2)  identified  secondary  phases, 

(3)  duplicate  detections,  (4)  regional  or  local  events,  (5)  velocity  failures, 

(6)  weak  signals,  and  (7)  data  dropouts.  Categories  (1)  and  (2)  contain  con-  . 

firmed  signals  and  the  rest  are  EP  false  alarms.  Definitions  of  these  seven 
categories  are  given  in  the  following; 
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(1)  Identified  events;  The  event  has  been  examined  by  the  analyst  and  con- 
firmed as  a P phase  of  an  event.  Events  In  this  category  are  published  In 
the  beam  location  summary. 

(2)  Identified  phases;  The  signal  has  been  confirmed  as  an  arrival  of  a 
secondary  phase  of  an  Identified  event.  Since  the  SDAC  bulletin  does  not 
report  events  without  first  confirming  P phase,  the  signals  In  this  category 
are  always  associated  with  category  (1). 

(3)  Duplicate  detections;  The  signal  Is  apparently  detected  by  a neighboring 
beam  (side  lobe  detections),  or  the  signal  Is  a coda  detection.  When  waveforms 
of  a duplicate  detection  are  displayed,  they  are  easily  recognized  and  rejected 
by  poor  signal  alignment  throughout  subarrays  or  they  are  obviously  part  of 
the  coda  of  an  event. 

(4)  Regional  or  local  events;  When  the  signal  characteristics  vary  distinctly 
from  subarray  to  subarray,  it  Indicates  that  the  signal  is  a regional  or  local 
seismic  event  arriving  at  LASA.  The  variation  in  signal  characteristics. is 
mainly  due  to  the  heterogeneity  of  the  crustal  structures  beneath  the  array, 

so  that  signal  coherencies  are  very  poor.  In  such  cases  the  signal  alignment, 
and  the  subsequent  attempt  to  define  the  beam  parameters  either  by  machine 
correlation  or  analyst  adjusted  alignments  are  not  reliable.  As  a result, 
the  analyst  rejects  the  event. 

(5)  Velocity  failures;  The  apparent  velocity  of  the  signal  is  higher  than 
the  theoretical  limit  of  P phase  velocity.  The  signal  could  well  be  a good 
signal  transmitted  through  the  Earth's  interior  core  from  a distant  location. 
However,  since  EP  is  not  presently  designed  to  recognize  core  phases  unless 
it  is  being  controlled  by  an  analyst,  this  type  of  signal  is  rejected  in 
automatic  processing.  Note  that  although  this  signal  is  rejected  because  of 
operational  restrictions.  It  can  be  very  useful  in  association  with  P phase 
detections  from  other  seismic  stations. 

(6)  Weak  signals;  The  signal  is  so  weak  in  amplitude  or  coherency  that 
neither  the  computer  nor  the  analyst  can  find  an  adequate  solution  to  define 
a beam.  However,  it  is  possible  for  an  experienced  analyst  to  recognize  the 
difference  between  a weak  signal  and  a signal  from  a local  event. 
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(7)  Data  dropouts;  The  detection  is  caused  by  bad  samples  in  the  data  stream 
(glitches)  or  the  detection  is  triggered  on  a sudden  data  dropout-restart 
situation. 

The  result  of  analyzing  all  EP  signals  during  this  test  period  is  summa- 
rized in  Table  I.  Identified  P phases  constitute  33.4%  of  the  totals  which  is 
comparable  to  the  result  of  an  earlier  study  made  by  Dean  (1972).  The  9.3% 
in  Category  2 is  the  result  of  the  analyst's  effort  to  identify  later  phases 
after  the  P phase  is  confirmed.  The  sum  of  these  two  categories,  42.7%,  are 
signals  from  confirmed  events.  The  remaining  57.3%  are  EP  false  alarms  of 
which  36.4%  are  due  to  duplicate  detections  and  11.6%  are  regional- local 
events.  There  are  no  false  alarms  due  to  noise  detections,  but  data  dropouts 
occurred  22  times  which  amounts  to  an  insignificant  0.8%  of  the  total. 

It  is  clear  from  Table  I that  duplicate  detections  and  regional- local 
events  are  the  dominant  causes  of  EP  false  alarms.  We  ask  whether  we  can 
reduce  them  by  simply  raising  EP  thresholds.  In  Figure  1 we  present  again 
the  seven  categories  of  EP  signals  Incrementally  grouped  in  signal- to-noise 
ratios  ranging  from  14  to  38  dB.  Percentages  of  each  signal  category  in  dB 
increments  are  computed  and  they  are  shown  in  a form  of  histogram.  This 
regrouping  shows  that  the  percentage  of  confirmed  events  steadily  increases 
as  the  threshold  is  raised.  However,  in  order  to  obtain  better  than  50% 
chance  of  confirmed  events,  the  threshold  must  be  raised  to  about  24  dB. 

The  analysis  of  Figure  1 shows  that  the  distribution  of  regional-local 
events  is  fairly  constant  throughout  all  S/N  ranges,  indicating  that  close 
range  events  are  frequently  detected  whether  or  not  there  is  a beam  directed 
toward  them.  The  rate  of  occurrence  in  each  S/N  range  is  approximately  10%. 
This  demonstration  clearly  Indicates  chat  raising  the  EP  threshold  will  not 
reduce  the  false  alarms  due  to  regional -local  events. 

Our  analysis  shows  that  there  are  pronounced  diurnal  variations  of  false 
alarm  rates.  To  demonstrate  and  Investigate  diurnal  variations  of  EP  false 
alarms,  we  regrouped  all  signals  according  to  local  time  of  the  day  at  LASA. 

In  Figure  2,  cumulative  hourly  frequencies  of  each  signal  category  are  tabu- 
lated in  terms  of  local  time  and  ovesented  in  two  histograms.  This  analysis 
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TABLE  I 


Classification  of  EP  Events  from  24  June  to  30  July  1974 

(Beam  Set  LBS*133) 


Category 

Number  of 
Events 

% 

(1) 

Identified  events 

958 

33.4 

(2) 

Identified  secondary  phases 

268 

9.3 

(3) 

Duplicate  detections 

1043 

36.4 

(4) 

Regional  or  local  events 

332 

11.6 

(5) 

Velocity  failures 

176 

6.1 

(6) 

Weak  signals 

69 

2.4 

(7) 

Data  dropouts 

22 

0.8 

Total 

2868 

100.0 

*LBS  = LASA  Beam  Set 


LEGEND 


Figure  1.  Percentage  distribution  of  EP  signals  in  increments  of  S/N 
ratio  for  the  period  24  June  to  30  July  1974,  when  LBS133  was  operative. 


LOCAL  TIME  OF  THE  DAY 


Figure  2a.  Diurnal  distribution  of  confirmed  signals  for  LBS133 
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showed  that  both  duplicate  detections  and  regional- local  events  are  highly 
concentrated  in  the  local  daytime  period  from  0800  to  1600  hours.  The  close 
correlation  of  both  categories  indicates  both  are  due  to  close  range  events, 
giving  us  the  first  concrete  evidence  of  the  relation  between  these  two 
categories  of  false  alarms.  In  order  to  reduce  the  number  of  these  false 
alarms  it  will  be  necessary  to  develop  some  method  for  monitoring  nearby 
events. 

The  remaining  question  is  the  cause  of  the  high  concentration  of  close 
range  events  in  the  local  daytime.  From  our  experience  in  operating  LASA  we 
know  local  mining  activities  can  generate  signals  that  result  in  detections. 
However,  whether  the  local  seismic  areas  such  as  Yellowstone  Park  area  or 
coastal  areas  of  Oregon-Callfornia  are  seismically  more  active  during  the 
daytime  is  not  known.  Since  the  signal  detection  is  based  on  S/N  ratio,  more 
signals  are  detected  during  the  quiet  nighttime  than  noisy  daytime  (Chang  and 
Seggelke,  1975) . This  effect  ia  demonstrated  in  the  diurnal  distribution  of 
confirmed  events  shown  in  Figure  2a.  In  this  figure  the  rate  of  confirmed 
events  is  higher  during  local  nighttime,  which  is  in  agreement  with  the 
result  of  Chang  and  Seggelke  that  showed  an  excellent  correlation  of  the  rate 
of  confirmed  events  with  the  hourly  noise  level.  We  do  not  assume  that  the 
overall  seismicity  within  the  surveillance  range  of  LASA  shows  diurnal  varia- 
tions as  was  suggested  by  Shlmshoni  (1971)  and  was  criticized  by  Flinn  et  al. 
(1972). 

What  remains  to  be  clarified  is  whether  particular  local  areas  do  or  do 
not  have  diurnal  seismicity  changes  that  can  be  related  to  the  sharp  increases 
in  both  regional-local  events  and  duplicate  detections  during  daytime.  In 
Figure  3a  we  have  plotted  the  number  of  daily  occurrences  of  duplicate  detec- 
tions and  regional- local  events  during  0800  and  1600  hours.  This  figure  shows 


Chang,  A.  C.  and  Seggelke,  R.  M. , 1975,  The  effect  of  band  pass  filters  on  LASA 
detection  performance;  SDAC-TR-75-9 , Teledyne  Geotech,  Alexandria,  Virginia. 

Shlmshoni,  M.,  1971,  Evidence  for  higher  seismic  activity  during  the  night; 
Geophys.  J.  R.  Astr,  Soc.,  v.  24,  p.  97-99. 

Flinn,  E,  A.,  R.  R.  Blandford,  and  H.  Mack,  1972,  Comments  on  ''Evidence  for 
higher  seismic  activity  during  the  night"  by  Michael  Shlmshoni;  Geophy. 

J.  R.  Astr.  Soc.,  V,  28,  p.  308-309. 


-16- 


that  these  two  categories  are  low  during  weekends  and  holidays,  suggesting 
that  local  disturbances  are  the  result  of  man-made  activities,  and  are  not 
due  to  diurnal  seismicity  variations. 

It  is  also  likely  that  codas  from  local  man-made  activities  will  raise 
the  ambient  noise  level  and  somewhat  impair  the  detection  capabilities  of 
teleseismic  signals.  In  Figure  3b  we  have  plotted  the  number  of  dally 
occurrences  of  identified  events  and  identified  later  phases  during  0800  to 
1600  hours.  A good  inverse  correlation  can  be  found  between  Figures  3a  and  3b 
in  that  the  number  of  identified  events  is  higher  on  Saturdays  and  Sundays. 
This  result  coincides  with  the  work  of  Woolson  (1976)  which  shows  that  LASA’s 
detection  threshold  (90%  confidence  level)  is  approximately  0.15  m^  better 
on  Sundays  when  compared  to  the  same  threshold  for  weekdays.  This  analysis 
clearly  shows  that  local  cultural  activity  is  a major  problem  in  the  LASA's 
detection  performance!  it  raises  the  rate  of  false  alarms  and  lowers  its 
detection  capabilities. 


Woolson,  J.,  1976,  LASA  detection  threshold  for  1974;  Comparison  of  Monday 
through  Saturday  with  Sunday:  Internal  Memorandum,  Teledyne  Geotech, 

March  1976. 
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NUMBER  or  EVENTS  NUMBER  OF  EVENTS 


Figure  3b.  Number  of  identified  events  (squares)  and  Identified  later 
phases  (triangles)  at  LASA  during  daytime  (0800-1600  local  time)  from 
24  June  to  30  July,  1974. 


DISCUSSION 


The  good  correlation  between  the  daily  occurrences  of  duplicate  events 
and  regional-local  events  shown  in  Figure  3a  indicates  that  duplicate  events 
result  mostly  from  regional-local  disturbances.  A good  reason  for  this  is 
that  there  are  no  regional-local  beams  in  this  beam  set.  Because  regional- 
local  events  are  generally  associated  with  large  amplitudes,  misaligned  signals 
are  detected  on  teleseismic  beams  as  side  lobe  detections  and  coda  detections. 

Because  the  real  interest  in  detecting  and  locating  events  with  large 
arrays  is  in  the  teleseismic  range,  and  also  because  of  operational  difficulties 
associated  with  detecting  and  locating  close  range  events,  past  operation  with 
LASA  has  not  used  beams  aimed  at  local  areas.  Regional  and  local  events  arrive 
at  LASA  with  relatively  low  apparent  phase  velocity,  and  many  beams  wJ.ll  be 
required  to  maintain  adequate  surveillance.  Even  if  the  signal  were  correctly 
detected,  because  of  its  high  frequencies  and  changing  signal  characteristics 
from  subarray  to  subarray,  it  would  be  difficult  to  locate  the  event  correctly. 
The  mission  of  the  array,  computational  time  and  core  requirements,  and  geo- 
physical difficulties,  discourage  attempts  to  treat  regional-local  events. 

These  are  the  reasons  why  the  DP  beam  set,  LBS133,  does  not  have  any  beams 
within  23  degrees  of  LASA. 

However,  the  results  of  our  analysis  in  the  previous  section  showed  that 
one  way  to  reduce  the  false  alarm  rate  is  to  eliminate  local  and  near  regional 
events.  Without  beams  in  close  range  areas  these  signals  are  detected  on 
teleseismic  beams  and  cause  and  Increase  in  the  number  of  false  detections. 

The  fact  that  we  cinnot  escape  detecting  local  events  leads  us  to  argue  that 
perhaps  local  areas  need  to  be  monitored  by  several  DP  beams.  By  detecting 
signals  on  close-in  beams,  side  lobe  detections  in  teleseismic  beams  can  be 
Identified  and  discarded.  Since  locating  nearby  events  is  difficult,  a simple 
algorithm  may  be  devised  to  eliminate  them  from  further  processing  by  limiting 
EP  to  teleseismic  signals. 

Adding  near  distance  beams  to  the  existing  beam  set  may  exceed  the  com- 
putational and  core  limits  in  the  computer,  and  this  aspect  must  be  discussed 
before  the  implementation.  The  old  concept  applied  to  the  LBS133  was  to  deploy 
beams  to  known  seismic  areas  with  equal  beam  separation  in  a hexagonal  pattern 
so  that  any  signal  from  a seismic  area  will  be  detected  in  one  of  these  beams. 
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The  consequence  of  rigidly  adhering  to  a pattern  is  that  few  of  these  beams 
are  exactly  placed  on  the  known  seismic  areas,  thus  all  DP  detections  are 
made  with  initial  location  errors.  If  DP  beams  are  selected  on  the  basis  of 
world  seismicity  and  are  aimed  directly  at  these  areas,  the  required  number 
of  DP  beams  can  be  reduced  and  the  initial  location  accuracy  can  be  increased. 
To  avoid  missing  events  from  non-seismic  areas,  coarsely  spaced  beams  can  be 
applied  as  a safety  precaution.  The  saving  in  the  required  number  of  beams 
can  be  used  for  local  beams  to  reduce  the  number  of  false  alarms. 

In  gereral,  side  lobe  detections  will  occur  when  the  peak  half-cycle  on 
one  instrument  is  added  to  the  following  or  preceding  peak  on  another  instru- 
ment in  the  beamforming  process.  This  type  of  false  alarm  will  occur  within 
one  or  two  cycles  of  the  main  peak  of  the  signal.  Since  teleseismic  signals 
are  known  to  have  a dominant  frequency  near  1 Hz,  side  lobe  detections  occur 
within  one  or  two  seconds  of  the  main  detection.  A detection  algorithm  with 
spatial  and  temporal  constraints  of  three  consecutive  threshold  crossings 
will  eliminate  this  type  of  false  alarm.  However,  during  the  analysis  we 
found  another  kind  of  large  scale  side  lobe  detection  associated  with  large 
signal  arrivals.  Suppose  as  an  example  a large  signal  has  just  arrived  at 
the  northern  most  subarray.  At  this  Instant  this  data  is  being  used  to  form 
beams  aimed  toward  the  south;  thus  a large  signal  in  one  subarray,  even  when 
reduced  by  a factor  of  N because  it  is  not  present  on  the  remaining  subarrays, 
will  produce  a false  detection  and  may  be  reported  as  an  event  arriving 
from  the  south.  Since  this  false  detection  is  spatially  and  temporarily 
separated  from  the  main  detection,  it  may  appear  as  a near  simultaneous 
arrival  of  two  Independent  events  from  opposite  directions.  Similar  detections 
con  be  obseirved  after  a large  signal  has  already  passed  through  the  array 
except  for  one  or  two  subarrays  at  the  edges.  The  result  of  a large  signal 
arrival  from  the  north  is  a set  of  three  detections  in  south-north-south 
beams  with  a few  seconds  separation.  For  an  array  with  50  km  in  diameter 
and  a large  signal  with  an  apparent  velocity  of  15  km/sec  =3.3  seconds 
before  and  after  the  main  detection.  Since  side  lobes  are  much  smaller  than 
the  main  detection  this  spurious  detection  pattern  can  occur  only  with  a 
large  signal.  We  think  that  the  pattern  of  detections,  the  size  of  the  main 
detection,  and  time  constraints  can  be  programmed  to  eliminate  this  kind  of 
false  alarm. 
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EVALUATION  OF  THE  NEW  BEAM  SET 


In  consideration  of  the  reasons  discussed  in  the  previous  section,  a new 
beam  set  was  designed  and  tested  with  the  on-line  DP  for  the  data  from  the 
period  of  June  13  to  June  30,  1975.  In  this  section  we  analyzed  and  compared 
the  performance  of  this  new  beam  set,  LBS160,  with  the  old  beam  set  LBS133. 

The  new  beam  set  has  a total  of  183  teleseismic  beams  aimed  at  locations 
for  which  we  have  LASA  travel  time  corrections  (Chiburis  and  Ahner,  1973). 

Since  these  locations  are  calibrated  for  travel  time  corrections,  we  expect 
these  beams  will  have  minimal  signal  losses  due  to  travel  time  errors.  In 
addition,  60  close  distance  beams  and  14  high  velocity  beams  were  selected  on 
the  basis  of  regional  and  PKP-range  seismicity.  For  the  purpose  of  covering 
non-seismic  areas  and  detecting  rare  events,  84  beams  are  added  to  the  beam 
set.  As  a result  LBS160  has  a total  of  341  beams  compared  to  300  beams  in  the 
LBS133.  Figure  4 shows  the  distribution  of  LBS160  beam  locations.  Detailed 
beam  parameters  are  given  in  Appendix  A. 

All  EP  signals  which  were  processed  during  the  period  of  LBS160  opera- 
tion were  analyzed  and  compared  in  the  same  way  as  the  previous  analysis. 

The  summary  of  signal  classifications  are  given  in  Table  II.  Comparing  the 
respective  categories  in  Table  I,  we  find  that  the  percentage  of  the 
Identified  events  has  increased  from  33%  to  42%  of  the  total,  and  Duplicate 
detections  have  decreased  from  36%  to  30%.  There  is  an  increase  from 
11.6%  to  13.4%  in  the  regional-local  events  presumably  because  the  new 
beam  set  has  a better  detection  capability  in  close  range  events. 

Next  we  conduct  a recount  of  all  EP  signals  shown  in  Table  II  with  a 
restriction  of  10  kra/sec  to  reject  all  regional-local  detections.  This 
restriction  eliminates  all  signals  detected  within  20  degrees  from  LASA,  The 
result  showed  21  identified  events,  161  duplicate  detections,  91  regional-local 
events  and  15  weak  events  were  eliminated  from  processing.  Table  III  shows 


Chiburis,  E.  F.  and  R,  0.  Ahner,  1973,  LASA  regional  travel  time  corrections 

and  associated  nodes:  SDAC-TR-73-6,  Teledyne  Geotech,  Alexandria,  Virginia. 
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Ld  map  showing  locations  of  LBS160  beams. 


TABLE  II 


Classification  of  EP  Events  from  June  13  to  June  30,  1975 

(Beam  Set  LBS*160) 


Catego 


Number  of 
Events 


TABLE  III 

Classification  of  EP  Events  from  June  13  to  June  30,  1975 
with  10  km/sec  Velocity  Restriction 


CateRorv 

Number  of 
Events 

% 

(1) 

Identified  events 

578 

51.5 

(2) 

Identified  secondary  phases 

87 

7.8 

(3) 

Duplicate  detections 

261 

23.2 

(4) 

Regional  or  local  events 

98 

8.7 

(5) 

Velocity  failures 

55 

4.9 

(6) 

Weak  signals 

15 

1.3 

(7) 

Data  dropouts 

29 

2.6 

Total 

1123 

100.0 
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the  percentage  distribution  of  each  signax  category.  Comparing  respective 
categories  in  Table  II,  we  find  the  percentage  of  identified  events  showed 
an  increase  from  42.4%  to  51.5%.  There  are  also  marked  reductions  in  false 
alarms;  duplicate  detection  from  29.9%  to  23.2%,  regional-local  events  from 
13.4%  to  8.7%,  and  weak  events  2.1%  to  1.3%.  The  result  proves  that  the 
velocity  restriction  is  Indeed  a very  effective  criterion  to  reduce  false 
alarms. 


Although  this  velocity  restriction  had  eliminated  a total  of  267  false 
alarms,  it  had  also  eliminated  21  good  local  events.  Among  these  21  events, 
we  found  one  signal  detected  at  16.5  dB,  three  detections  between  18  to  20 
dB,  and  17  remaining  events  with  better  than  20  dB  S/N  detection.  It  is 
therefore  possible  to  set  a higher  threshold  at  about  18  dB  to  process  local 
detections  occurring  within  20  degrees  from  LASA. 


By  excluding  local  events  there  are  a total  of  578  teleseismic  events 
during  the  17  day  period,  or  an  average  of  34  events  per  day.  This  is  the 
highest  daily  average  we  have  ever  obtained  at  SDAC,  In  Figure  5,  recurrence 
curves  for  LBS160  are  shown  together  with  similar  curves  for  LBS133  computed 
for  the  year  1973.  There  are  a total  of  8197  events  during  334  days  of  1973, 
or  a dally  average  of  24.5  events  per  day.  Detection  thresholds  for  these 
two  recurrence  curves  were  computed  with  two  methods:  the  first  with  the 

maximum  likelihood  method  by  Ringdal  (1975),  and  another  by  fitting  a best 
linear  line  estimate  between  magnitude  ranging  from  3.8  to  4.9  nij^.  By  using 
the  maximum  likelihood  method,  90%  detection  thresholds  does  not  change  much, 
nij^  » 3.88  for  LBS133  and  m^^  = 3.84  for  LBS160.  However,  the  90%  detection 
thresholds  are  quite  different  if  we  fit  a straight  line  only  to  the  upper 
portion  of  the  recurrence  curve.  We  find  nij^  = 3.8  for  LBS133  and  m^^  = 3.55 
for  LBS160.  It  can  be  seen  in  Figure  5 that  the  maximum  likelihood  method 
has  a better  fit  over  wide  range  of  magnitudes,  but  the  linear  line  estimate 
is  better  for  the  specified  range.  We  think  the  true  detection  threshold 
change  lies  somewhere  between  these  two  values.  Since  local  events  are 


Ringdal,  F.,  1975,  On  the  estimation  of  seismic  detection  thresholds:  Bull. 

Seism,  Soc,  Am.,  v.  65,  p.  1631-1642. 
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excluded  in  this  comparison,  we  conclude  that  there  is  some  indication  that  the 
more  accurate  and  v;ell  calibrated  teleseismic  beams  have  improved  detection 
capability  by  the  order  of  0.1  + .1  m^^. 

Figure  6 shows  the  performance  of  LBS160  in  terms  of  percentage  distri- 
bution of  EP  signals  grouped  in  the  seven  categories  specified  earlier  as  a 
function  of  incremental  S/N  threshold.  The  most  significant  improvement  in 
this  figure  is  the  disappearance  of  regional-local  events  in  high  S/N  ratio 
detections.  Tliis  means  that  although  difficulties  with  small  signals  still 
exist,  most  of  large  local  events  can  be  properly  processed  and  identified. 

The  accuracy  of  the  teleseismic  beams  can  be  evaluated  by  comparing  travel 
time  errros  of  the  DP  beams.  The  travel  time  errors  are  differences  of  the 
final  travel  times  of  identified  events  (i.e.,  the  travel  time  associated 
with  the  final  location)  and  travel  times  of  the  DP  beam  which  detected  the 
signal.  Smaller  travel  time  errors  indicate  DP  beams  are  more  accurate  and 
thus  associating  DP  detections  with  another  station  will  be  optimized.  Figure  7 
shows  the  comparison  of  travel  time  error  is  -14.385  seconds  for  LBS133  and 
-3.733  for  LBS160.  Standard  deviations  are  28.56  for  LBS133  and  22.96  for 
LBS160.  An  obvious  skewness  can  be  observed  in  the  distribution  curve  of 
LBS133,  but  no  obvious  skew  can  be  seen  for  LBS160.  We  believe  LBS160  is 
generally  superior  to  LBS133. 

In  summary  there  are  three  reasons  that  implementation  of  new  beam  set 
and  false  alarm  criteria  will  lead  to  improve  performance  at  LASA.  The  first 
reason  is  the  reduction  of  travel  time  errors  in  the  new  beam  set-  The  mean 
travel  time  error  for  the  old  beam  set  is  ~14  seconds,  which  is  equivalent 
to  an  average  initial  location  error  of  200  kilometers  at  A = 65  degrees. 

The  mean  travel  time  error  for  the  new  beam  set  is  ~4  seconds  and  that  is 
equivalent  to  an  average  initial  location  error  of  50  kilometers  at  the  same 
distance.  Secondly,  new  beam  locations  are  well  calibrated  for  travel  time 
anomalies  and  may  well  therefore  have  a lower  detection  threshold.  Thirdly, 
the  reduction  of  false  alarms  can  reduce  the  computer- analyst  workload  so  that 
we  can  lower  the  operating  threshold  without  difficulties. 
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CONCLUSIONS 


As  a result  of  investigating  on-line  seismic  signal  detections,  we 

found  that  seismic  signals  from  local  and  regional  events  cause  most  false 
alarms,  and  seismic  noise  causes  few  false  alarms.  Codas  from  these  close 
range  events  also  disturb  the  detection  capability  of  the  array  by  masking 
small  teleseismic  signals. 

Local  and  near  regional  events  are  technically  difficult  to  confirm; 
however,  we  found  that  steering  beams  to  close  range  areas  somewhat  reduces 
the  rate  of  false  alarms.  In  order  to  discriminate  against  local  events, 
detection  beams  must  be  deployed  into  these  areas.  This  arrangement  will 
reduce  side  lobe  detections  and  coda  detections  on  teleseismic  beams. 

Although  the  deployment  of  local  beams  alone  can  reduce  false  alarms, 
we  found  the  most  effective  criterion  is  to  set  a velocity  restriction  to 
eliminate  these  signals  detected  on  local  beams.  Since  most  local  events 
are  detected  with  high  S/N  ratios,  an  arrangement  to  set  a high  S/H  threshold 
on  local  beams  will  pass  good  events  and  eliminate  false  alarms. 

Instead  of  using  several  beams  equally  spaced  in  a particular  seismic 
area,  we  found  it  more  effective  to  place  one  beam  directly  on  the  center  of 
the  area.  Comparison  of  recurrence  curves  showed  that  such  a new  beam  set 
has  a lower  detection  threshold.  This  is  perhaps  due  to  two  factors;  detec- 
tion beams  are  closer  to  true  epicenters  (thus  less  beamforming  loss),  and 
these  beams  are  well  calibrated  for  travel  time  anomalies  for  the  particular 
areas  used.  Deployment  of  new  seismic  beams  showed  that  the  average  location 
errors  of  DP  detections  are  reduced  from  200  km  to  50  km. 
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